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ABSTRACT 


The conventional low-power CW doppler projectile tracking radar 1s investigated 
through simulation. The power requirement for full range tracking with real-time track 
data processing stopped within 160 m of the point of impact is obtained. This power 


can be lowered bv a factor of more than 4 if a set of three switchable filters is utlized. 
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I. INTRODUCTION 


At a proving ground, it is always desirable to record the trajectories of test-fired 
munitions for later analvsis. This is especially true When anomalies occur. Continuous 
wave (CW) radars have often been adopted for this purpose because of their simplicity. 
Compared to a pulsed radar of sinnlar range performance, a CW radar has no modulator 
in its transmitter. hence no range gating 1s necessary In its receiver. Furthermore. a CW 
radar has a bandwidth corresponding to the doppler shift, which rarely exceeds 100 
kilo-Hlertz (KHz) while a pulsed radar has a band width in the mega-Hertz range. Thus 
a CW radar can be designed to carry out coherent signal integration over a longer time 
interval. This increases the integration gain and lowers the required transmitter power 
for the CW radar. 

High gain antennas of verv narrow beamwidths are also used to reduce the trans- 
mitter power requirement. It 1s important to aim such antennas accurately on the target. 
To acquire the bearing and elevation angles of a target directly, several receiving feeds 
together with complicated comparator and multiple data channels are needed. Since the 
trajectory of a test round can be roughly predicted, an alternative to the actual acquisi- 
tion of target tracking angles is by steering the antenna according to the predicted tra- 
jectory. The directly measurable data will be the doppler shift of the target as a function 
of flight time. Discrepancies between this data and the doppler shift based on the pre- 
dicted trajectory should be carefully analyzed and explained. When an anomaly occurs 
which causes the target to break off track, this recorded data will be valuable for diag- 
nosing the probleni. 

Many CW radar tracking svstems are available. For example, a model DR 582 X- 
band Instrumentation Radar manufactured bv TERMA Elektronik AS of Denmark 1s 
being installed at the Ta-Fu Proving Ground in Taiwan. These systems typically track 
a long-range projectile only up to the summit of the trajectorv. Since the point of inipact 
of a long-range projectile is of great interest, it is desirable to extend the range of the 
radar to cover as large a portion of the trajectorv as possible while keeping the svstem 
relatively simple. The purposes of this thesis are to investigate the maXimum range at- 
tainable using a simple CW radar and to determine the required minimum transnutter 


power. 


When tracking a projectile, a CW radar detects only the doppler frequency shift due 
to the projectile motion. This frequency shift is proportional to the component of the 
projectile velocity pointing radially away from the radar. By combining this information 
with the tracking angles, the ground range of the projectile can be deduced through in- 
tegrating this velocity component. A low-pass filter can be employed to perform this 
signal integration in real time. The output can be used to update the predicted track if 
a data filtering mechanism such as the Kalman filter is utilized in real time. Since in- 
formation on the doppler integrating filter is not provided bv any manufacturer, one way 
to implement this low pass filter 1s suggested in Chapter 2. The design described therein 
1s consistent with the operational instruction for the aforementioned system. It is needed 
for tracking simulation required in Chapter 1V. The trajectorv-updating Kalman filter 
is a desirable improvement to currently available svstems on the market. Its design 1s 
outside the scope of this thesis. 

As explained above, trajectory estimation 1s necessary for tracking with a simple CW 
radar. Since the radar and the gun are moving with the spinning Earth. the familar 
reference frame in which the radar and the gun are stationary 1s an accelerating frame 
Which 1s difficult to work with. Instead, by ignoring the rotation of the Earth around 
the sun, the inertial frame of reference having its origin at the center of the Earth and 
Its Z-axis along the spin axis of the Earth simplifies the equations of motion of the 
projectile. The z-x plane can be chosen to contain the radar at t= 0 when the projectile 
is fired. To obtain the familar ground-based quantities. transformation to the inertial 
frame in Which the radar (or the gun) 1s passing instantaneously through the z-x plane 
can be performed. In this frame, the distance and the relative velocity of the projectile 
from the radar (or the gun) can be found easily. 

The projectile will be spinning around its axis if a rifled gun barrel is used. The lift 
force due to air pressure will rotate the projectile to align its axis along its air velocity 
vector. The kinematics of the projectile, including its angular momentum, ts most easily 
described in terms of the motion of, and the rotation around, its center of mass [Ref. 
1]. The orientation of the projectile, which lies along one of its principal axis of inertia. 
can be described 1n terms of two of the three Euler angles [Ref. 2], leaving the third angle 
for measuring the spin of the projectile. 

It can be shown that the effect of gravitational force on the projectile orientation 1s 
negligible. The ratio between the gravitational force on the center of mass to the force 
affecting the angular momentum is of the order of the square of the ratio of projectile 


size to the Earth radius [Ref. 3]. Thus the change of the orientation of the projectile axis 


to 


is due largely to the resistance of air. For a streamlined body. the hft from the aur flow 
is large when there is an angle of attack {Ref. 4]. Furthermore. the rotation of the 
projectile axis is slow. It 1s reasonable to assume that the orientation of the axis of the 
projectile is always along its air velocity. This assumption will be utilized to determine 
the radar cross-section of the projectile. 

It is easily seen in the ground coordinate that the projectile velocity 1s much greater 
than the wind velocity so that the latter only introduces small corrections to the 
projectile axis. Since the wind acts on the projectile over the entire trajectory, a better 
estimate of the range can be obtained by including the wind velocity profile to trajectory 
prediction. A less elaborate alternative is to adopt a wind profile model such as the 
Ekman spiral for wind speed approximation [Ref. 5]. Since the change in the predicted 
range is not expected to be large, the wind velocity 1s treated as constant throughout this 


thesis. 


II. PROPOSED CW RADAR DESIGN 


A. RADAR SYSTEM SCHEMATICS 
The receiver design of a radar is critical to its range performance. Since no specific 


information for any CW radar tracking system is available, the following receiver con- 


figuration 1s proposed. Here Figure 1 shows the schematics of the radar. 
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Figure 1. Schematics for Proposed CW Radar Design. 


The received signal at the frequency (f,+,;) from the receiving antenna 1s fed to a 
narrow-band, low-noise RF pre-amplifier with a bandwidth of 1 MHz centered at 10.475 
GHz. The noise figure of such a pre-amplifier can be as low as 1.5 to 2 dB, with a 
munimum gain of 28 dB [ Ref. 6 ]. This low-noise pre-amplifier can limit the receiver 
noise figure ( F, ) to about 2 dB. The | MHz bandwidth 1s chosen according to the es- 
timated maximum doppler frequency ( f,....) and the fact that the target 1s alwavs flving 
awav from the radar. The output signal from the low-noise pre-amplifier is reduced to 
the base band at the mixer, and then passed through the video filter of bandwidth &,. 
The filter output, at the doppler frequency f{, can be recorded either with an analog or 
a digital recorder. The bandwidth of this recorder should be wider than /;,,, to allow 
truthful recording of the raw data for later analvsis. The doppler signal is then fed into 
the square-law detector. 

The output of the square-law detector oscillates at a frequency of 2f, Since this 
frequency is proportional to the velocity of the projectile receding away from the radar, 
it can be integrated to obtain the projectile range information. 

The integration is most easily accomplished by counting the number of cross-over 
points of a pre-determined level by the square-law detector output. As to be explained 
later, this number, denoted by MN, is the displacement of the projectile awav from the 
radar measured in half the transmitter wavelengths. 

The square-law detector, a counter, and a clock will be adequate for measuring the 
range and the velocity of the projectile receding from the radar. First, a measurement 
base, MB, is set in terms of the number MN of transmitter half wavelengths. Then the 


time 7, for the output of the square-low detector to run through MN pre-set level 
MIB 
if 


nf 





crossings 1s measured. Here MB is the radial range increase of the projectile and 
is the radial velocity of the projectile. 

As discussed in the introduction, the trajectory of the projectile is predicted before- 
hand. The elevation and bearing angles of the radar are updated at the end of every 7, 


interval. This information is sent to the servo control system to keep the antenna on the 


mroyectile. 


B. THEORY 

The theory of operation for the proposed CW radar system will be described. The 
doppler effect is used in a CW radar for moving target detection. This effect was found 
by Christian Doppler in 1842. It states that the frequency of a harmonic signal will shift 


when reflected by a moving target. The change in frequency is proportional to the ve- 


locity of the target. Let f, denote the doppler shift, , the radial velocity of the target 
with respect to the radar. c the velocitv of hght, and f the transmitted carrier frequency. 


The doppler frequency shift is 


V, 

gee (2.1) 
Si 

= 21,54 (2.2) 


The doppler is zero when the trajectory of the projectile is perpendicular to the radar 


line of sight. It 1s positive for an approaching target (see Figure 2). 
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Figure 2. Spectra of Received Signals: (a) no doppler shift, no relative motion (b) 


approaching target (c) receding target. 


The frequency of the signal out of the square-law detector 1s 2/,. Thus the number 


of crossings, MeN. of this signal over a pre-set level Within @ time Intemiain] ace 


“fy | 
ae (2.3) 
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According to equation (2.1), the radial displacement of the projectile over this period 
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Figure 3. NMfeasurement Base and Velocity vs. Time. 


MB = T,xV. (a 


fs 
“wee” 


= MVx+ (2.6) 
Pemexdinpie, 1 the transmitting frequency /, of imterest is 10525 MHz, we can set 


the number MN to 100. Then 


x 
mp = 100% 2.907 x 107 (2.7) 


Peloses < 10: 
= 1.4234m (2.8) 


This is a small change in target location for updating the radar pointing direction. 


Since iB is the displacement of the target along the radar pointing direction, this dis- 





placement can be projected onto the ground. The ground range cal 
ing up the individual projected ranges. 


II. PROJECTILE TRAJECTORY PREDICTION 


To analyze the range performance of the radar, information about the trajectory of 
the projectile is required. There are many factors such as air drag. gravitational force, 
projectile shape and spin, Earth rotation,...etc., which can affect the trajectory. As dis- 
cussed in the INTRODUCTION, the projectile spin mainly provides a stabilizing 
mechanism to the projectile and 1s assumed to have little effect on its trajectory. The 
direction of the projectile axis (or heading) is assumed to be along the direction of the 
air Velocity of the projectile. Thus only the motion of the center of mass of the projectile 
meeas to be considered. The projectile center of mass and the center of air pressure are 
assumed to be the same. The trajectory is thus reduced to that of a particle trajectory. 

Wind velocity can increase, decrease or deviate the trajectory through air drag on 
the projectile. In this thesis, it is assumed that wind velocity 1s constant over the course 
of the projectile; it 1s usually called the ballistic wind. 

Drag force depends on the air velocity of the projectile. It points against the 
projectile heading. Its magnitude is proportional to the square of projectile air speed. 
the local air density. the projectile cross-section, and the drag coefficient. Air densitv 
@epends On the temperature profile of the atmosphere. The drag coefficient is verv 
complicated and, in fact, varies with the projectile air speed. Measured data of the drag 


coefficient are available. The drag force on the in-flight projectile is given below : 
rs pau 1 yay? 2 
= — A> ACD — Hh) erl) 


where H¥" is wind velocity, I’ is the velocity of the projectile, C, 1s the drag coefficient, 
A 1s the cross-section area of the projectile. and p is the local air density. Here h is the 
heading of the projectile. which is a unit vector in the direction of aay 

The gravitational force is the only other factor affecting the trajectorv in an inertial 
Heme Of reference. In the accelerating reference frame fixed on the surface of the Earth. 
there is the fictitious Coriolis force. It causes an apparent trajectory drift. As observed 


on the ground, this drift can reach 3.] km for a trajectory of 27.4 km ground range. 


A. THE TRAJECTORY IN A FIXED INERTIAL FRAME 

In order to simplify the equation of motion, the Earth-centered inertial frame (r, 
8,@) is used. The plane containing the radar and the spin axis of the Earth when the 
projectile is fired at t = QO 1is the z-x plane. The radar is inoving with the Earth at an 
angular velocity of Q. Its location in this inertial frame is (7,. Op, 922) . 


The components of the projectile velocity in this reference frame are : 


r=! (328 
rd=J% (33a) 
r(sin @)d = Vs (3.4) 


The Lagrangian formulation of the principle of least action is convenient for deriv- 
ing the equations of motion of the projectile (Ref. 7]. The Lagrangian of the projectile 
in the inertial spherical coordinate system [Ref. | , 8] (r, @. ), without considering the 


drag force, can be expressed as : 


L = T(r, 6. 6,7, 0.6) — U(r, 6, 6) (3.5) 
Vil el, fel eee ee aaa 
L= Tl POG PO S112) a aera) (3.6) 


where ¢ is the gravitational acceleration at ground level. For the r component, the 


equation of motion is: 





GE 

sili, se a 
Gaec? Cr oe 
: a ee 
Vee CO a se (3.8) 

r 

For the 6 component : 
a. EE CL 
ds Se 3.9 
at =a CO ot 
@ (726) = 76? sin @ cos 0 3.10 
Vp) = en COs (3.10) 


and, for the @ component: 


d ¢L Gl 

Se ee aa EE ————— a 

ai db CH oe 
ua ai”  sin’é) = 0 (3.12) 
dt a ms 


The nonlinear drag force can be included. The complete equations of motion become: 


Nps 

== eee. p 

a —(-)g+—— +h, (3.13) 
me 7) = tel tof 3a 
pea EON e * Ne (3.14) 
i. | aaa 
aU sin @V,) = rsin 6f. (als) 
r=acth (3.106) 


where h 1s the height from the Earth surface and ais the radius of the Earth. The ratio 


of projectile height to the Earth radius can be expressed as : 


r h 
a as (3.17) 


aq 


Miesprojcctile velocity is : 


ee 1 +104 1',0 (3.18) 
= prt +00 +ysin bad (sly) 


and the wind velocitv 1s : 
WV = W641 64107 (3.20) 


Note that /,, f,. f, are the components of the drag force divided by the mass of the 


projectile : 





ee, ——— AC. | VV = it" |? (3.21) 


LAT 


=f t+ foot feo (3.22) 


The equations of motion can be discretized. Assume that the gravitational and the 
drag force acceleration are approximately constant over the period 7, when the projectile 


moves from r, to r,_,, the equations of motion are approximately : 


Ponmaa? 2 
n+] n Qa 2 ae 2 3 
7 = = Sener +r, sin 0,0, + fen (3.239 
a n 

2 e ? e 

Pie ar 

n+) nt] n ae Me ; - 
a = sind, cos 6,62 + raho, (3.24) 


T,, 


ee ; Zs 
Fe esi Gein, — Fein as 


7 


Qn 


= 7, 8in Of, (3.25) 


From the above equations, if v,. 6,. @, and r,. @,. @, are known, then 7, Gea 


M4. U4. and Gecam pe den\ceasiollonus: 


With measured /_. equatronn( ly wrewides je 4 then 


6.4, = 6,+— >—-T, (3.2% 
Equations (3.27) and (3.23) give @,.,. then 


d, a ee 
ee = P, a | Vhs (3.257) 


Note that the muzzle location is r, and the muzzle velocity of the projectile, which 
is assumed to be known. is J’. Through the above procedures, the full trajectory of the 


projectile can be constructed bv computer simulation. 


B. GROUND-OBSERVED TRAJECTORY 
In order to obtain the trajectory as observed from the ground, the position of the 
projectile should be converted into its relative location from the radar antenna and from 


the gun. The necessary coordinate transformation matrices are given below (Figure 4 


on page 13). 





Figure 4. Inertial Coordinate Svstems. 


1. Trajectory Relative to Radar 
To set the elevation and bearing angles for the radar and to estimate the doppler 
shift, the location and velocity of the projectile relative to the radar as viewed from the 
ground are needed. On the ground, the projectile motion is most easily described in an 
inertial frame in which, at time, t, the radar 1s passing through the z’ — x’ plane instan- 
taneously. The 2’ axis is chosen so that it coincides with the z axis. The transformation 
from the fixed inertial spherical coordinate system to this (x’, »’ ,z’) svstem can be ac- 


and.) : 


complished through the following rotation matrices R,, R a 


rec } 


COS OT Wisin Oe 


R(Qr) = | —sinQ2r cosQr 0 (3228) 
0 0 ] 
cos Pp QO  —<sin P, 
Krel@p) = Si @p 0 COs, (3.30) 
0 | 0 


Sep co ae 
Rene) —aecOs.0,, sii Grae (3.39) 
0 0) | 


a. Velocity Transformation 


Ul 


The projectile velocity m the (x’, 4”. z’) coordinate svstem can be obtained 
by first changing the results mn section II]A from the spherical coordinate system into 


the rectancular coordmate systema ye 


V px Le 
Voy | = Rreldp)RekOn)} Moe (3.32) 
I ps I po 


Since the (x’, 3”, 2’) coordinate svstem 1s obtained by rotating the (jaya 


coordinate svstem by an angle {7 around the z-anis. 


Sse Ls 
oe = RAO ae (3.33) 
ve, Ve 
ie 
= Riedy ~ QDR, A8,)| Uno (3.34) 
ie 


The velocity of the radar, | ,, in thel(x 73 42 )coordmaresrr 


Vy! = Rye (9) Ry AER) 0 (3.35) 


sin@p cos@p Of; 0 


= 0 0 | 0 (3.36) 
cos@p, —sinOr Of Pp, 
0 
= er (oro) 
0 
Pence 
Ve = Vr" (3.38) 
= (rp sin Op)" (3:39) 
b. Projectile Location Relative to the Radar 
iiiterelationship between (x. 3 .2 and (1, 0, @) 15 given by : 
x rsin@cos @ 
aie— | r 5in@ sind (3.40) 
z HiCOS @ 
x’ x 
Vo] = RLQn! 3 (3.41) 


i 
hy 


rsin@cos(@ — Qr) 
=] rsin@sin(d — Q2) (3.42) 


RiCOSO 


c. Radar Coordinate 
Substitute (7,, @,. 821) inte equation Se) 


foes rp Sin Op 
Vp == 0 (3.43) 
Ir rp COs Op 


2. Trajectory Relative to Gun 
Relative to the gun (muzzle), a similar inertial frame (x ~ .y ~~ ,z im Whiehagae 


gun 1S instantaneously passing through the z”-.x” plane at time t can be used. The 


z” axis coincides with the z-axis. The coordinate transformation is : 


So rsin @ cos(@ — dg — 21) 
ve" | = [ rsin@ sin(d — dg — Q (3.44) 
ri rcosé 


The firing plane of the gun is the plane containing its barrel and passing through 
the center of the Earth. Vor artillery. it 1s customary to use the coordinate svstem 
(x..4,.2,) In which the gun muzzle is the origin and the z,— x, plane 1s the firing plane: 
Here 2,18 the vertical direction and x, 1s along the horizontal from the gun. The relation 
between (7,3 2,) and (x ~ .3'" ,2z “) depends on the gun bearing angle @,)5. winichwisaaae 
angle of the x, axis measured counterclockwise from the East. If the the gun is oriented 


along 0,4 and is located at r, . @¢. @,. the vector r, — 7, im the coordinate (x, 3,2) oe 


X/ O. “COS:0Ues sy sinigye. cos@g 0 sn@g ffx” 
yp = | OP si Ga eeCOS Dy aa ex 0 I 0 io (3.45) 
s I 0 0 —sin@e 0 cos@g |} z” 


Note that +, — 7, is the relative location of the projectile from the gun. Hence 


its components can replace (x~ ,;° ,2 +) directly injequatiens ey 


16 


IV. RANGE ANALYSIS AND POWER REQUIREMENT 


A. RANGE EQUATION 
The power requirement of the CW radar can be obtained from the radar equation 
met. 9] ; 


Sf) 


CoE Cree are Sta 
P= 2 (4.1) 
G or 
In equation (4.1), the radar cross-section (o). the maximum range (R), and the 
propagation factor (F) will be discussed in this chapter in detail. The Boltzman constant 
femeees 7.0) n-my ke- A. The absolute temperature J, 1s 290°A . The wavelength / 1s 
0.0285m. Other parameters are based on the configuration of the CW radar as proposed 
Mee napter 2. Thev are defined as follows. 
1. The Receiver Noise Figure 
The receiver noise figure (/,) 1s the ratio of the signal-to-noise power ratio at the 
receiver input to the signal-to-noise power ratio at the receiver output, assunung that the 
feeetver temperature is 290°A . The F, value for the receiver considered is assumed to 
2. System Loss Factor 
The loss factor consists of the fractional loss of transmitter power delivered to 
the transmitting antenna (L,). the fractional loss of the receiving antenna output power 
Selmeered to the receiver input (/,). and miscellaneous loss (L,). They are discussed be- 
HE SN7. 
a. Transmission Loss 
The transmission loss is the loss in the transnutting or the receiving portion 
of the transmission line. Since the transmitter 1s light-weight and 1s to be mounted verv 
close to the antenna, L, 1s estimated to be about 1.2 dB due to the use of narrow band- 
width rotary joints. The receiver pre-amplifier can be mounted directly to the receiving 
antenna output, followed with the mixer immediately. This arrangement will reduce L, 


to about the same level as L.. 


b. Miscellaneous Loss 
The miscellaneous loss (L,) for this CW radar system is mainly due to the 
quantization loss at the comparator. Since only two-level (crossing) comparison is 
made, it 1s not expected to be significant and will be ignored. 
Base on the above argument. the total loss is estimated at 2.4 dB. 
3. Signal-to-Noise Ratio 

For a tracking radar, it 1s important that the receiver output signal is sufficiently 
high above the noise level so that information carried by the signal can be extracted re- 
ably. With the proposed receiver configuration. the signal-to-noise ratio (S N) of the 
output from the video filter is the figure of significance. First of all, this signal is re- 
corded for later analvsis. Its minimum S N has to be assured. In real-time. after passing 
through the square-law detector, this signal generates the timing sequence for updating 
the radar pointing direction and for evaluating the projectile radial velocity. It drives the 
gain control to keep the reference signal at a preset level. It is also the input into the 
comparator following which the number of signal cross-overs is tallied and the elapsed 
time 1s clocked as explained in section IIA. 

The square-law detector and the comparator perform threshold detection. The 
reference signal acts as the threshold. The comparator outputs a +1 signal when the 
detector output is above the reference level: it outputs a -J signal when the detector 
output falls below the reference. The reference signal level, which 1s set at a fractionsem 
the detector output signal, determines the accuracy of the real-time output of the svstem. 
In what follows, the minmum S N rauo and the threshold setting in terms of the 
threshold-to-signal power ratio are discussed. 

The signal input to the comparator is oscillating at twice the doppler frequency. 
Within each cvcle of the signal, the portion rising above the threshold level to the peak 
and back to the threshold level can be considered as a pulse. This pulse 1s contaminated 
with noise. which 1s assumed to be white. The probability of each contaminated pulse 
falling under a certain threshold level is well-known. This probability of missed de- 
tection 1s the fraction of signal pulses failed to be detected by the comparator. Between 
two pulses, the number of false alarms 1s the number of extra cvcles counted. These two 
effects tend to cancel each other so the optimum threshold setting 1s achieved by equal- 
izing the fraction of missed pulse and the fraction of noise pulses. The program to 
compute these fractions are listed in Appendix D [Ref. 10]. 

Figure 5 on page 19, Figure 6 on page 20, Figure 7 on page 20, Figure 8 on 


page 21, and Figure 9 on page 2! show the fractional errors due to missed pulses and 
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noise pulses respectively. If the signal-to-noise ratio is 10 dB, the comparator output 
error 1s 5.55% when the threshold-to-signal power ratio is set at 0.455. Stronger signal 
allows the threshold to be set at a lower fractional value with an improved error ratio. 
For example, for a signal 12 dB above the noise, the threshold setting of 0.38. or -4.2 
dB below the signal, will limit the error to 1.15%. In this thesis. the minimum signal- 


to-noise ratio 1s chosen to be 10 dB. 
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Figure 5. Missed Detection and Noise Error vs. Threshold-to-Signal Power Ratio 
S N= 10 GB. 
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Figure 6. Mflissed Detection and Noise Error vs. Threshold to Signal Power Ratio 
SS = 105 dh. 
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Figure 7. Missed Detection and Noise Error ys. Threshold to Signal Power Ratio 
S N= I! dB. 
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Figure 8. Missed Detection and Noise Error vs. Threshold-to-Signal Power Ratio 
S N= 11.5 dB. 


Missed Detection (+) & Noise Error (0) 


0.42 0.44 0.46 


Threshold-to-Signal Power Ratio 





Figure 9. Missed Detection and Noise Error vs. Threshold-to-Signal Power Ratio 
S N= 12 dB. 
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4. Transmitting and Receiving Antenna Gain 
Separate but identical transmitting and receiving antennas are used to reduce 
direct coupling of the transmitter power into the receiver. Since this separation is small, 
it is neglected in the analvsis. In x-band, antennas with one-degree beamwidth are easilv 


achievable. Based on this beamwidth, the antenna gain (G) [Ref. 9] 1s : 


ae 
eg (4.2) 
AQ 
= 44.3dB (4.3) 
and the diameter D of such an antenna 1s 
D = 72x (4.4) 
6 
— 
aoe — (4.5) 


S. Video Bandwidth and Filters 

The required minimum video bandwidth (B,) is a function of the doppler fre- 
quency. As previous mentioned, the bandwidth should be greater than /,,,, . Although 
the doppler frequency f, of the projectile increases rapidly mght after firing. 1t decreases 
monotonically afterwards bv a factor of more than 4 over its trajectory. An example 1s 


shown in Figure [0 on page 23. 
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Figure 10. Doppler Frequency vs. Slant Range. 


Referring to the range equation (4.1), in order to use the lowest transmitted 
power possible for full range tracking. the radar range and video bandwidth can be 
traded off. Note that. assuming a constant radar cross-section, the transmitter power 
depends on the product R*B.. A filter of narrower bandwidth mav be emploved when 
the projectile reaches a greater distance. Adoption of a set of three switchable doppler 
filters 1s recommended. The following is an example for its implementation based on the 
doppler frequency of the projectile given in Figure 10. 

The maximum doppler shift /,,. from the projectile shown in Figure 10 1s a 
little less than 54 KHz. The maximum range R is about 27 km. Three filters are to be 
used over three different ranges. Each of the filters will have a bandwidth of 1.2 times 
the maximum expected doppler frequency over the range in which it 1s to be used to al- 
low for a 20%0 estimation error. The transmitter power level will be minimized if the 
fourth power of the maximum range times the maximum doppler shift in the corre- 
sponding range is a constant. 

a. First Filter 

When the projectile 1s within 16.06 km slant range, the video bandwidth 1s 
maeisen tO be 1.2 x f...,. Thus: 
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B, — | Xi (4.7) 
= 63.9 kH: (4.8) 


b. Second Filter 
When the projectile moves bevond 16.06 km. but within 23.73 km slant 
range, the doppler frequency has a maximum value of about 22.03 kHz. The video 
bandwidth ts : 


B= 0a (4.9) 
= 26.44 4H: (4.10) 


c. Third Filter 
When the projectile reaches 23.73 km slant range, the maximum GOppieg 


frequency 1s about 12.7 kHz. The video bandwidth is : 
Be 2 ol ee (4.11) 
= 1552 hee (47125 


Since A*f,.,, 1s now less than I 4 of the original value when onlv one video filter 1s used, 


the munimum required transmitter power will be reduced by more than a factor of 4. 
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B. TRAJECTORY SIMULATION AND RESULTS 

Based on the procedures for trajectory estimation presented in Chapter III, a simu- 
lation program Was written and 1s included as Appendix A. It was written in PC-Matlab 
386. The program is described as follows. 

1. Input and Parameters. 

Gee aril 
The operation is assumed to take place in Taiwan. The radar 1s Jocated at 

the tongitude of 121°36’27 ” East and the latitude of 24°36'38 ” North. The Earth ra- 
Gus iS approximately 6.3677225x= 10m (average value of the artic radius 


6.356912 x 106 m and the equator radius 6.378533 x 106 m). The gravitational acceler- 


ealome at sea level is 9.81 m/sec’. The angular velocity of Earth spin 1s 
fee 21 1S44 =x 10~ rad sec [Ref. I 1}. 
b. Radar 


The transmitting frequency is 10525 MEiz. The measurement base MB can 
be varied to change the update rate of the antenna pointing direction. For both the re- 
ceiving and the transnutting antenna, the beamwidth 1s 1° and the gain is 44.35 dB. The 
antenna height and the distance between the gun and the radar can affect the iniual 
readings of the doppler signal. 

c. Gun 

The firing elevation angle is assumed to be 45° and the azimuth 1s counter- 
clockwise from the East and 1s also assumed to be 45°. Both the elevation angle and the 
azimuth angle can be varied. The distance between the gun and the radar can be ad- 
justed. (Variations of up to 2000 m are possible at the test site.) The barrel length of 
the gun under consideration 1s 7 m. 

ad. Projectile 

The target is a 155 mm cahbre artillery projectile. Its weight 1s 45.4 kg, with 

a total length of 938 mm. The maximum muzzle velocity is 850 m sec. 
e. Environment 

(l;) Temperature and air density profiles. The temperature 7, at sea sur- 

face 1s assumed to be 288°A and the temperature at the projectile height h is based on 


the following profile [Ref. 12]: 


ee 1, — ay (4.13) 


where h is the height from sea level, a; is the temperature lapse rate (-0.0065 °K m). 
The temperature (7,) is used for calculating the air density (p,) at the projectile: 


Niele aan ass 
Ph = ie) Gere (4.14) 
O 


In equation (4.14) p, 1s the air density at sea level (1.225 kg'm?), and R jis the gas con- 
stant (287.05 n-m kg-°A). 

(2, Wind profile. The wind is assumed to be constant during the flight 
time over the projectile trajectory. The value used is the so-called ballistic wind. It can 


De ex presseqias ; 
W= WF ++ Wd (4.15) 


In this thesis. the wind is 6 m'sec, North-West. It should be noted that the extent to 
Which the estimated trajectory can be modified by the actual wind profile is subject to 
further study. 

(3, Mach number. Vhe Mach number 1s the ratio of the projectile air 


speed to local sound speed. The sound speed at the height of the projectile 1s : 





a, = . ART, (4.16) 
Where Kk 1s 1.4 for idea gas. The Mach number (M) is given by : 


i Us a (4.18) 
Ch 
where | 1’ — Wy is the projectile speed with respect to air. 
f. Drag Coefficient 
The drag coefficient (C,) 1s required for calculating the drag force) Itwdes 
pends on the shape of the projectile. The C, curve is measured experimicnta! it eeeee 
given as a function of Mach number M. For a 155mm projectile. the C, curve is siiamam 
in Pisuve 7) somepdce: 27. 
2. Results of Simulation 
The results of trajectory simulation are described below. The variable names 


used in the program are given in the parentheses when they are first menuoned. 
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a. Radar Observed Trajectory 
The slant range (SRR) is the distance from the radar antenna to the in-flight 
projectile. The curve showing the projectile height (Height) versus its slant range from 
the radar 1s shown in Figure 12 on page 28. The total flight time 1s 87.4 seconds, and 
the maximum range 1s 27.4 km. 
b. Gun Based Trajectory 
The ground range 1s measured from the gun to the projectile. The projectile 
height vs. in-flight time (t) and the projectile height vs. ground range are shown in 
Figure 13 on page 28 and Figure 14 on page 29. The maximum ground range (Grange) 
is 27.4 km. 
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Figure 11. Drag Coefficient vs. Nlach Number for 155 mm Projectile. 
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Figure 12. Projectile Trajectory: Height vs. Slant Range. 
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Figure 13. Projectile Trajectory: Height vs. Time. 


28 


F 
= 
wee” 
eet 
sox 
ame 
© 
a 
v 
Oo 
— 
c. 


1S: 
ground rance (m) 





Figure 14. Projectile Trajectory: Height vs. Ground Range. 


The cross range (CRange) measures the drift of the projectile from the firing 
plane. The firing plane is the vertical plane containing the gun barrel. In the case con- 
sidered, the maximum cross range 1s about 6.3 km South-East of the firing plane. The 
point of impact 1s about 3.14 km away from the firing plane. The curve of cross range 


vs. ground range 1s shown in Figure I5 on page 30. 
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Figure 15. Projectile Trajectory: Cross Range vs. Ground Range. 


c. Radar Tracking Parameters 
The azimuth (traz) and elevation (trel) antenna pointing angles are required 
for tracking. The curves of these angles vs. in-flight trme when the gun-radar distance 
is 500 m are shown in Figure 16 on page 31 and Figure 17 on page 31. 


The azimuth and elevation angular speeds of the antenna are shown in 


Figure J$ on page 32, and Pieure 19 70n pagers”. 
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Figure 16. Radar Azimuthal Tracking Angle vs. Time. 
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Figure 17. Radar Elevation Angle vs. Time. 
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Figure 18. Radar Azimuthal Tracking Angular Speed vs. Time. 
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Figure 19. Radar Elevation Tracking Angular Speed vs. Time. 


The antenna servo system has to point the antenna to within half a beamwidth in the 
direction of the projectile. Thus the half beamwidth divided by the angular speed of the 
antenna is the required response time of the servo system. The servo response time Vs. 
the projectile in-flight time 1s shown in Figure 20 on page 33. From the above results, 
it can be seen that over the first few seconds, the elevation angular speed is verv high. 


The distance between the gun and the radar can be increased to reduce this speed. 
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Figure 20. Servo System Response Time vs. Time. 


d. Doppler Frequency and Velocity 
The doppler frequency (f;} is a directly measurable data by a CW radar. 
The doppler cvcle can be detected with the square law detector. In the proposed CW 
radar configuration of this thesis, the doppler cvcle is derived from the time interval T, 
for each measurement base traveled. This timing sequence has to be simulated in the 


’ 


) is calculated 


dopp/ 


tracking program. In the simulation program, the doppler velocity (J 
from the projectile velocity relative to the radar and the radar pointing direction. Then 
T, 18 estimated from the time required to travel the distance MB with this doppler ve- 
locitv. The curves of doppler frequency vs. slant range, doppler frequency vs. time, and 
doppler velocity vs. time are shown in Figure 10 on page 23, Figure 22 on page 34, and 


Figure 2] on page 34. 


a8 


Based on the above results from trajectory simulation, the projectile radar 


cross section and the transmitter power can be estimated. Thev will be discussed later. 
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Figure 21. Doppler Frequency vs. Time. 
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Figure 22. Doppler Velocity vs. Time. 


C. PROJECTILE RADAR CROSS-SECTION 
The projectile radar cross-section (RCS) is denoted by o. It is the area intercepting 
the amount of plane wave power which, when scattered equally in all directions, 


produces an echo at the radar. In other words, the RCS ts given by : 
5 
a Anr?( = - (4.18) 


@eeere ©. 1s the Scattered electric field strength from the target. r is the distance from the 
target in its far field where the quantity £, is measured, and E, is the incident plane wave 
field strength at the target. 

Determination of the RCS of a complex target is difficult. Usually the radar cross- 
Section of a target depends on the following factors: a) target size, b) target shape, c) 
radar frequency. d) target aspect angle, and e) radar polarization. The RCS 1s greatly 
simplified if the target size is large compared to the radar wavelength. Consider the ra- 
dar cross-section of a simple sphere. The RCS as a function of the circumference of the 
sphere measured in wavelengths (2za/z) is shown in Figure 23 on page 36. The region 
Sameve the tarcet size 1s larse compared to the radar wavelength, 1.e., where 2za// 1s 
greater than 10, 1s called the optical region [Ref. 9]. The RCS in this region 1s roughly 
the physical cross-section at area of the target. 

fire 155 mim projectile is in the optical region at the transmitter frequency 
(10525 MHz) of interest. Its RCS value will be approximated with its geometrical 
cross-secuonal area. Because of the rotational svmmetrv of the projectile, the radar 
cross-section of the projectile is a function of the angle @,., between the radar pointing 
direction (D) and the projectile heading (h). 


The radar pointing direction D is CeneDy 
(4.19) 


Where rp is the projectile position vector, rz is the radar position vector as defined in 
Sirapter 11]. 
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Figure 23. Radar Cross-Section of the Sphere vs. as : 





The heading of the projectile, h, is in the direction of its air velocity. It 1s also de- 
fined in ChapterlII: 


In terms of h and D, the aspect angle @,., of the projectile as seen from the radar 1s 
Gasp = cos '(D . hr) (4.29 


Note that 0 < @,,.< x. For simplicity, the projectile is considered as a circular cone on 
top of a circular cylinder, as shown in Figure 24 on page 37. Tl.e height of the cylinder 
is 414.4 mm. The height of the cone is 523.6 mm. The diameter of the cylinder and the 


base of the cone 1s 155 mm. 
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Figure 24. Approximate Shape of 155mm Projectile. 


Seae vertex angle of the cone, 2a, is given by: 


me x aE G Ore 
= (2.1469 rad. (423) 
= 8.42° (4.24) 


The phvsical cross-section of the cone is the projected area of the projectile along D. the 
radar pointing direction. onto a plane perpendicular to D. This cross-sectional area will 
be discussed below according to whether | Deh —=2€0S SGT | Deh = COS.00 
l. Case | 
As | Deh |=cosoa, the projection of the projectile along D is shown in 


feieure 25 on page 38. 


of 





Figure 25. Radar Cross Section Case 1. 


The total RCS value 1s 
o,= na? |D-h|+2ab| Dxh| (4.25) 


where ain the radius of the cvlinder and 6 is the height of the cylinder. 
Die Heo aSee 
As | Deh |< cosa, the total RCS is given by equation (4.20) (Figure 26)6n 
page 39). 


projection plane 
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Figure 26. Radar Cross Section Case 2. 
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(4.26) 


(4.27) 


Co A re (4.31) 


A 4 A 
A, DOS 0 — J) COs. Ore 
SS eee (4.32) 


oe, 4 ee 2 
~ Sin G55, COS! O +> SI) Gage yee 


where L is the height of the cone and NV is the normal vector of the surface A,. The de- 
tailed derivation of the RCS equations is given in Appendix B. 

From the above results. the RCS value of the projectile in-flight can be pre- 
dicted. It 1s shown in Figure 27 and Figure 28 on page 41. From above figures, it can 
be seen that the maximum RCS 1s about 0.09 >”? and the minimum RCS is about 
0.02 27. Since the aspect angle of the projectile changes verv fast right after it is fired, 


the RCS value varies rapidly at the beginning of the trajectorv. 


“ 
t 
& 
~~ 
© 
£ 
© 
b. 
« 
2 
“ 
c 
c 
oe 
u 
t 
A 
" 
" 
© 
& 
u 
b. 
€ 
~ 
L -} 
x 


36 468 


Aspect angle (degree) 





Figure 27. The In-Fligkt Projectile Radar Cross Section vs. Aspect Angle. 
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Figure 28. In-Flight Projectile Radar Cross Section vs. Slant Range. 
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D. PROPAGATION EFFECTS 

Propagation of radio waves through the Earth’s atmosphere is significantly affected 
by the propagation medium. The presence of the Earth complicates the situation in 
three wavs: (a) the surface of Earth scatters radiation and produces an interference pat- 
tern, (b) the Earth casts a shadow and gives rise to diffraction phenomena. and (c) the 
Earth’s atmosphere is inhomogeneous and produces refraction effects [Ref. 13]. It is 
customary to group all these factors into a single quantity called the propagation factor 
(F). In this thesis, the Engineer’s Refractive Effects Prediction system (ERE?S)siiem 
]4] is utilized to compute the propagation factor based on local environment. 

EREPS uses rav tracing for radio Wave propagation within the horizon. It includes 
the bending effect of the refractive index variation of air and the divergence factor when 
the ray 1s reflected from the Earth. The field at the target 1s the sum of the directa 
and the reflected ray. resulting in interference peaks and nulls which become significant 


when the target falls near the lienizom(Tieures. |; 
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Ficure 29. Reflection and Interference from a Spherical Earth. 


For the over-the-horizon region, the Waveguide mode approach is adopted by 
EPs Phe Earth § atmosphere is modeled as a leaky waveguide. The Earth's surface 
is the lower boundary of the guide. Modes of high attenuation along the horizontal di- 
rection are called leakv modes. Modes of low attenuation are called trapped modes. In 
x-band. the effect of the evaporation duct is significant. Trapped modes occur when the 
evaporation duct height is large or at least comparable to the radar wavelength (sce 
Figure 30 on page 44) [Ref. 15]. 

1. Computation 

A typical EREPS output is given in Figure 30 on page 44. The local annual 
average meteorological conditions for propagation are available {Ref 13]. The following 


parameters are used : 


J. Evaporation duct height eli) Joe! 
fk (efieemue Earth radius factor) Sey) 

3. Wind speed wo iS 
4. Absolute humidity 3 WO ieee 
Ee ouriace refractivity ee? 


Note that the projectile is at a particular height when it arrives at a range. This 
computation has to be repeated for each point along a projectile trajectory. This 1s done 
for several trajectories at various ranges (Figure 32 on page 45 and Appendix C ). Also 
mete that, in Figure 32 on page 45, tle ground range is given a negative sign because it 
is measured from the point of impact toward the radar. 

2. Summary of Results 
a. Signal Strength Enhancement 
PAOMMEICUre 52 Om pace 45, it can be Seen that, with ducting, which is 
present most of the time at the test site, the return signal will be substantially enhanced 
from the corresponding free space value. (Note that the computation 1s carried out for 
a target at | m height. 1.e., when the projectile is near its point of impact. The applicable 
mae fOr the projectile of interest lies between 27 to 28 km.) Thus ducting is helpful in 
enhancing the quality of the recorded full range tracking data. 
b. Strong Interference 
As shown in Figure 32 on page 45, the return signal will exhibit strong in- 


terference when the projectile reaches about 250 m within the point of impact. The more 
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than + 15 dB fluctuation makes it difficult to track the projectile in this region using the 
simple CW system proposed in this thesis. On the other hand, since this strongly fluc- 
tuating region is within a few hundred meters from the impact point, this small portion 
of the trajectory can be estimated with good accuracy. Hence it is desirable to select a 
‘stop-processing point’ at which the CW real-time radial velocity and range computa- 
tions will be halted. The radar antenna will be guided by the predicted trajectory, and 


the return video signal will still be recorded. 
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Figure 30. Waveguide Mode with Ducting Height and Antenna Height. 
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Figure 31. Propagation Factor with 13.1m Duct Height. 
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Figure 32. Propagation Factors vs Ground Range from Impact Point. 
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3. Stop-Processing Point 
It should be noted that EREPS does not take into account the antenna radi- 
ation pattern. Since the antenna to be used have a beamwidth of I° only, interference 
between the direct ray and the reflected rav will not be important until the antenna ele- 
vation angle is lowered to within 0.5° above the horizontal direction (Figure 33 on page 
46). 
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Figure 33. Narrow Beamwidth Effect on the Spherical Earth Surface. 


This criterion can be used to determine the stop-processing point as follows: when the 
system records an antenna elevation angle falling below 0.5°, the square-law detector 
and the comparator will be switched off. The servo will steer the antenna according to 


the predicted trajectory smoothly toward the projectile point of impact. Characteristics 
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for two such stop-processing points for a radar-to-gun distance of 500 m and two dif- 
ferent radar heights of 10 m and 30 m respectively are displayed in Table 1 on page 47. 
Note that these point are close to where the propagation factor starts to fluctuate 
pmronmeis as predicted bv ERPS. 

Since the processing of real-time data 1s stopped before contributions from the 
reflected ray become important; in what follows, the propagation factor will be consid- 


ered as unity. 


Table 1. STOP PROCESSING POINT. 


















Antenna height (m) 
Radar to gun distance(m) 
Projectile geieoe oa 


x 
Elevation angle at one point (degrec) 
ee veUlon angle difference (degree) 
Azimuth aatie at Rana point( degree) 
Azimuth angle diflerence to impact point(degrce) 
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FE. POWER REQUIREMENT 

From equation (4.1), the power requirement for tracking up to the stop-processing 
point is shown in Figure 34. In order to reduce the required transmitter power re- 
quirement, a multiple filter design is proposed for this system. 

The multuple filter design 1s based on the characteristics that, after an initial rise, the 
doppler frequency decreases with increasing slant range (Figure 35 on page 49). By 
switching to video filters of narrower bandwidths at longer ranges, the coherent inte- 
gration time of the radar is increased hence the power requirement can be lowered. 
Combining the fixed range independent parameters in the radar range equation. the 


power requirenient’ can bere pressed acs 





7 a R°B, 
P, = (3.7339 x 107**) x (—s— Day (4.33) 
Re 
= Cx (Sa ere aia (4.34) 
C = 4.4807 x 10°” (4.35) 
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Figure 34. Power Requirement for Full Range Tracking without Multiple Filters. 
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Since f,.¢ and R are functions of the trajectory (Figure 35 on page 49, and Figure 36 
on page 50). if the tracking range is divided into three regions 1, 2, and 3, and one filter 


is used for each region, the transmitter power requirement becomes : 


4 
P, = C( “GC )y max (fa) max (4.36) 
Ro 
= C( ml ae Ge oenas (4.37) 
ae 
= C( “GC Vice’ Grscae (4.38) 


where the constant C js the same as that of equation (4.34). 
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Figure 35. Doppler Frequency vs. Slant Range. 
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Figure 36. Radar Cross Section vs. Slant Range. 


The power requirement written in terms of that of the single filter case 1s: 


where C,. the power reduction factor, is the factor by which the power réequiremenmae 


reduced. 


In order to find the optimal switch-over points for the filters, the logarithm of f, 
versus K*/o, normalized to (/,),,,, and (R*/¢),., respectively, ts plotted (Figure 37 cippaae 
il): 

Assume that the first switch-over point is located at (x,, 4), and the second one ts 
located at (x,,3;). These values can be determined from Figure 37 on page 51 by ob- 


serving that, at the optimal switch points : 
eae (4.40) 


= (4.41) 
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The optimal locations so obtained are listed in Table 2 on page 52. The power re- 
duction factor is C, = 10°*!* = 4.186, as seen from Figure 37 on page 51. 

The power requirement as a function of slant range when three filters are used 1s 
plotted in Figure 38 on page 52. The antenna height is 10 m and the radar-to-gun dis- 
tance is 500 m. The required power is 37 Watts. A conventional CW tracking svstem 
will track to the summit point only. This point is located at a range of region | and re- 
quires a transmitter power of 19.5 Watts. Almost full-range tracking is possible by 


doubling the conventional power through the use of multiple filters. 
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Figure 37. Log(Normalized f;) vs. Log(Normalized R‘/c). 
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Figure 38. Power Requirement for Trajectory Tracking with Multiple Filters. 


Table 2. MULTIPLE FILTERS SWITCH POINTS 


Parameters region | region 2 region 3 stop track- 
ing point 
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fod 
ae 
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Slant range } Slant range (km) | — 


td 
= 
oe 


oa 
Maximum / (KHz) eee eS 


Video bandwidth 63.90 22 

(KHz) 

Radar cross section 0.0385 (0.0438 0.0862 0.0903 
()77" ) 

Required power 36.29 5625 36208 

(Watts) 
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V. CONCLUSIONS 


It is demonstrated in this thesis that full-range projectile tracking, with real-time 
track processing stopped at a range within 160 m of the projectile point of impact. can 
be achieved if the transmitter output power is 37 watts or higher and a set of three 
switchable filters are used. If a single filter is used, an output power of 160 watts 1s re- 
quired. Either configuration can be implemented with a commercially available, small 
size, air-cooled RF amplifier which can be mounted near the transmitting antenna. 

The tracking relies on the predicated trajectory. Although real-time doppler infor- 
mation 1s available from the radar, it 1s used only to confirm that the projectile 1s on 
track. A algorithm to utilize this real-time information for updating the predicted tra- 
jectory is desired. Kalman filtering appears to be a promising technique which deserves 


further investigation. 
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APPENDIX A. PROJECTILE TRAJECTORY SIMULATION AND 
RANGE ANALYSIS 


clear 

t2=clock; 

format long e; 

7 the first capital letter designations are: P for the projectile; 
7% G for the gun (muzzle); R for the radar. 


% The numerals 1, 2, 3 are for the radial, theta and phi components in 

% an inertial spherical coordinate system; 4 is for the radial component 
vo in the cylindrical coordinate system associated with the spherical 

% coordinate. The origin is at the center of the Earth; the polarwgz) 

% &@Xis is the axis of spin of the Earth. 


7% East and South designations are local quantities in a frame rotating 
vo With the Earth 


% constant declaration 


EarthT=86164. 06; %o Sidereal period of Earth spin 
EAV=2“pi/EarthtT; % Earth angular velocity 
A=(65569 127637050 3)7 2: % approximate Earth radius 
g=9. 81; % gravitational acceleration 
RLATI=24+36/60+38/3600; % latitude 24036'38" (radar site) 
RLONGI=121+52/60+27/3600; *% longitude 121052 27" (radar sagem 
f=10525E6; % radar transmitted frequency Hz 
lambda=3E58/f; % radar transmitted wavelength m 
“input parameters 
{del ww 
diary ww 
RB=input('radar to gun base distance(m) =? 
% RBaz=input(' radar to gun direction from east=? ')*pi/180 
RBaz=0; 
° MN=input('Measurement base number setting =? >) 
MN=1000; 
Rh=input( ‘Radar height(m)=? 
° Gel=input(' firing elevation angle(degree)=? ")/180*pi 
Gel=45/180*pi; 
% Gaz=input('firing azimuth angle(degree)=? ')/180*pi 
Gaz=45/180"pi; 
diary 


answer=input('If you want reenter parameter please enter 1 else enter ); 
if answer==1 


return 
XXa 
end 
B1=5; % gun barrel length 
A=pa- COL 572). % projectile cross section area 
m=45. 4; %o projectile mass (kg) 
FAC=0. 5*A/m; % coefficient for computing drag 


54 


MB=MN**( lambda/2); % measurement base 


% Radar Location at t=0 % the inertial z-x plane 
Rl=a+Rh; 
R2=(90-RLATI)*pi/180; 
R3=0; 
R2s=sin(R2); 
R2c=cos(R2); 


RG4=R1*R2s; 
RZ=R1*R2c; % z-coordinate of radar 
RV=R4“EAV; % radar velocity 


% Firing Direction and Speed 
PVo=850; % muzzle speed of projectile 
Gels=sin(Gel); 
Gelc=cos(Gel); 
Gazs=sin( Gaz); 
Gazc=cos( Gaz); 


% Gun (muzzle) Location at t=0 


Gh=B1l**Gels; % muzzle height 

Gl=a+Gh; 

G1=B1*Gelc; 

RGs=-RB*sin( RBaz)-G1l*Gazs; %® distance, south of radar 
RGe=RB**cos( RBaz)+G1*Gazc; % distance, east of radar 


G2=R2+RGs/G1; 
G4=G1*sin(G2); 
G3=R3+RGe/G4; 
Gz=G1*cos(G2); 


% Displacement Transformation 

G2c=cos(G2); 

G2Zs=sin(G2); 

Pee vont yoAZo, J, -GAZs ,GAZC; 1,0,0=*-G2C,0,G28; 0,1,0; -G28,0,G2C); 


yeecoyectile Initial Location 
P1=G61; 
Pin=Gl; 
P2n=G2; 
P3n=G3; 
P4n=G4,; 


% Projectile Initial Velocity 
PV1n=PVo*Gels; 


PVg=PVo*Gelc; %© Velocity along ground 
PVe=PVe*Gazc; % velocity, east 
PAV2n=-(PVg"Gazs)/G1; % angular velocity, theta 
PAV3n=PVe/G4+EAV; % angular velocity, phi 


% Drag Force Considerations 


% Wind Profile: assumed constant in this program. 
% The geostrophic wind should be measured and the Ekman spiral 
7% model can be adopted for a better approximation if the true 


La) 
" 





o/ 


7 Wind profile is not available. 


% Wind=input('Enter the wind speed m/sec at muzzle height’ ) 


Wind=6; 


o/ 
/o 


wind velovity at muzzle height 


* WindBearing=input(' Enter the true wind bearing in degrees’) 


WindBearing=135; 


vo 


Waz=((90-WindBearing)+180)*pi/180; 


We=Wind*cos(Waz); 
Ws=-wWind*sin(Waz); 


clockwise from neren 


% Wind velocity, east 


Temperature Profile (assumed linear) 


% wind velocity, souen 


% Th=input('Enter the height (m) at which temperature is measured’ ) 


Th=10; 


7 TREF=INPUT( ' ENTER THE MESURED TEMPERATURE IN CELSIUS’ ) 


Tref=15; 


Trefk= 
Thi=Th+ta; 
S =-0. 0065: 


% lP-TREPR So lin), % 


APALRK. DENSI 


% DH=INPUT( ENTER THE HEIGHT (M) AT 


DH=GH; 


Clos bonlrer. 


» TEMPERATURE LAPSE RATER Ti 
AIR TEMPERATURE AT THE PROJECTILE 


WHICH AIR DENSITY IS MEASURED’ ) 


% DAO=INPUT('ENTER THE MESURED AIR DENSITY IN KG/CUBIC METER’ ) 
% Air density at sea level 


Dao=l. 225; 


TDao=TrefKk+s**(Dh-Th); 


SS—=207.05¢ 


KC=1. 4; 
RAg=-1l-g/(s*ss); > 
% DA=DAO*(TP/TDAO )**RAG; ” 


7 opeed of Sound 
KR=KC*ss; 


% Vah=sqrt(KR*TP); 


”% gas constant N-m/kg-K 


*. Gib © Neto rare te 


o 


% Calculation faceou 
AIR DENSITY AT PROJECTILE HEIGHT 


% speed of sound at Pl 


% curve fitting for Drag coefficient(experimental data aaa 


x=( 

6507.6. 

es ee 52.9); 
y=(0. 16,0. 16,0. 16 
0: 16,0716. 0. 16.0. 
02875 ,0).200n0 27 


cds=spline(x,y); 


ree 


0,0.05,0.15,0. 2,0. 23,0. 25,0.3,0.35,0 4 (04550 ere 
Te. 0. 825, oF 875 AIG 98, Leis 05, ye eis 2, ee S i alee 


0.16,0.16,0.16,0.16,0. 16,0. 16,0516 (GC 


- 162,0. 17,0. 205,0. 225,0. 27,0. 282,0. 29, 
25,0, 23 0.195505 ie 


% Trajectory Computation Begins 


OY whe alone oie ala ale wlan’ ente ale ale aloats ole 
#u fe ah Fe Oe Ft FX HA FH OF C8 OF OD FH 


while Pin>a, n=n+tl; 


% Update 
Pi=- ine 
P2=P2n; 


erie? earl eri eris 


o/ 


main loop begin #*%vrssdrsede 


spline curve fitting funeerem 
get spline coefficient 


initial time 


sevesesdese secs sesesestesodesesevetedesecesk 


P3=P3n; 
P4=P4n; 
PV1=PVi1n; 
PAV2=PAV2n; 
PAV3=PAV3n; 
PV2=PAV2*P1; 
PV3=PAV3*P4; 
PVe=( PAV3-EAV)*P4; 
Height(n)=Pl-a; 
% velocity transformation matrix 
PR3=P3-EAV*t(n); % location relarive to radar 
PG3=PR3-G3; % location relative to gun 
P2s=sin(P2); 
P2c=cos(P2); 
PR3s=sin(PR3); 
PRoe=—cos( PR3); 
PG3s=sin( PG3); 
PG3c=cos(PG3); 
ViRPI=P2s*PR3c; 
VTR21=P2s*PR3s; 
VTR31=P2c; 
pee ey OR, P 2G PR3c, -PR3s; VIR2Z1 ,P2c*PR3s ,PR3c; VIR31,-P2s ,0); 


VTG11=P2s*PG3c; 
WIGZ1=P2s*PG3s; 
NEGsl=P2c; 
% Relative Location 
Perey 2—( PIl*=ViRI1; VIRZ]; VIR31)-(R4: 0; RZ); 


SLR=NORM(PRXYZ); ‘7o SLANT RANGE FROM RADAR 
eo oN 
Wes His hacho rhc isekerhtriee arise Wes eis hse hg tai ber ise hi eh te br arr kar ir) ark ha rh te tr) oo he oir de a Be) She) Se Soh 
DRC: e- =PRxyz/SLR: di eos relatuwe to radar 
e—OR( 1 n)*R2stDRC3,n)*R2c; % tracking direction, elevation 
Dn=DR(3,n)*R2s-DRC1,n)*R2c; PmeEceranecearrecrion, north comp. 
e=DR(2,n); % tracking direction, east component 
traz(n)=atan2(Dn,De)*180/pi; 7e azimuthal angle 
ee. Ae ee ee % elevation angle 
ote BE TEBE TE ETE EET ETE EE TEE TE ric chacigelichs clocks celta ylshaachae arse de lee hair] arial oie taieh dr) oe) Oe oie oe ie ae he ee eae hed 
PGXYZ= DTG (PL¥- VTG11: VTG21: WiGol)-(64-- GZ): 
SL=norm(PGxyz); % Slant range 


SRange(n)=SL; 
GRange(n)=a*acos((P1**2+G1**2-SL 2)/(2*P1*G1)),; 
% ground range 
PGx=PGxyz(1); 
PGz=PGxyz( 3); 
FRange(n)=a*atan( PGx/(G1+PGz)); % range in firing plane 
CRange(n)=PGxyz(2); % cross range from firing plane 


*Relative Velocity 
fe—(PV1, PV2; PV3); 
PVR=VTR*PV-(0; RV; 0): 
sn )=PVR; % velocity relative to radar 
Wine (nl )=PV 1; 


7% Time Interval and Doppler effect Estimation 
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% 


or 
4/0 


VDopp=PVR' *PRxyz/SLR; 
Vdop1(n)=VDopp; 
fd1(n)}=2*VDopp/ lambda; 
dtn=MB/VDopp; 

£( t= et) en: 
Frictional Acceleration 


TP=[Tretkts(Pl-tn >. 
Da=Dao*(TP/TDao)**RAg; 
Vah=sqrt(KR*TP); 
PVW=CPV1; PV2-WS; PVE-WE); 
PVWM=norm( PVW); 
PVWMM=PVWM/Vah; 
PH=PVW/PVWM; 

HeadingR(: ,n)=VTR*PH; 
CD=ppval(cds ,PVwhn); 
FA=-FAC*Da*CD*PVWM*PVW; 


Location and Velocity Prediction 


PA 1=(PV2"**2EVS 2) Flt EAC) —2 Ga ba) 


PAA2=P2s*P2c%*PAV3**2+FA(2)/P1; 
PAA3=FA(3)/P4; 

PVin=PA1*dtn+PV1; 

Pin=0, 5*CRVIFPVin)“dtn- er: 
PAV2n=( PAA2*dtn+PAV2)*(P1/Pin)**2 
P2n=0.-5*{ PAY 27 PAVen Gather 2. 
P4n=Pin*sin(P2n); 

PAV3n=( PAA3**dtn+PAV3)*(P4/P4n)**2 
P3n=0. 5*( PAV3+PAV3n )*“dtn+P3; 


end 


& 
/o 


° 
£0 


- a's a! os ala -" ° i 
a'se's& 
es raved ee ri eek eu ee 


Estimate Point of Impact 

Th=P1l-a; 

IA=-PAI]; 

IV=-PV1; 
dt=(sqrtClV**2+2*ih* lag Ty 7 TA 
N=nt+1; 

CON =t(n dz: 

Pin=a; 

PV1in=PA1*dt+PV1; 
PAV2n=(PAA2*dt+PAV2)*(P1/Pin)**2; 
P2n=0. 5*( PAV2+PAV2n )*dt+P2; 
Poen=Pingsint’ Zn). 

PAV 3n=( PAA3*dt+PAV3 )*(P4/P4n)**2; 
P3N=0. 5**( PAV3+PAV3N )*DT+P3; 


Update 
P1=Pi1n; 
PZ=P2n; 
P3=P3n; 
P4=P4n; 

PV 1I=PVin: 
PAV2=PAV2n; 
PAV3=PAV3n; 
PV2=PAV2*P1; 


e 
bd 


D) 


5§ 


te 
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Doppler velocity 


% Doppler Frequency 
square-law doppler detector 


air temperature at the projectile 
air density at projectile height 
speed of sound at Pl 

PROJECTILE AIK VELOC 


projectile heading 
heading as seen from radar 


drag acceleration 
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se end of main loop souk desea sk seor see scat it 


PV3=PAV3*P4; 
Height(N)=0; 


% velocity transformation matrix 


PR3=P3-EAV*t(N); ~) 1 ©edtion relative to radar 
Pe3-PR3-G3; % location relative to gun 


P2s=sin(P2); 

P2c=cos(P2); 

PR3s=sin(PR3); 

PR3c=cos(PR3); 

PG3s=sin(PG3); 

PG3c=cos(PG3); 

VTR11=P2s*PR3c; 

VTR21=P2s*PR3s; 

weno l=P2c; 

em Pee 2C--PR3C.=-PR3s3 VIR2Z1,P2C*PR3S,PR3C; VIR31,-P28,0); 


VTG11=P2s**PG3c; 
VTG21=P2s*PG3s; 
VTG31=P2c; 


memelative Location 
PRXYZ=(P1*-VIR11; VTR21; VTR31)-(€-R4; 0; RZ)-; 
pe norm ERxyzZ); 


y Ai ee Radar parameters GEGES Ema LOM Me<ecrerrr ects Zea ever res 2s 
DR(: N= =PRxyz/SRR(N); ‘0 location relative to radar 
B2=DR( 1, N)*R2s+DR(3 ,N)*R2c; wetcacking Girection, elevation 
DN=DR(3,N)*R2S-DR(1,N)*R2C: APIRACKING DIRECTION, NORIH COMP. 
De=DR(2,N); % tracking direction, east component 
traz(N)=atan2(Dn ,De)*180/pi; % azimuthal angle 
oN be Oe eon s ee Bee 
yet en ees ents ee 7 oO Or ht igri caches Sor hiriate btiie he hrc Goan nics Siro cckas 
PGxyz=DTC* (PL: (VTG11: VTG21: VTG31] -{ G4: 0: Gz]): 
SL=norm( PGxyz); % Slant range 
SRange(N)=SL; 
GRange(N)=a*acos((P1**2+G1**2-SL 2)/(2*P1*G1)); 
% ground range 
bex=—PGxyz( 1); 
PGz=PGxyz(3); 
FRange(N)=a*atan(PGx/(G1+PGz)); % ground range in firing plane 
CRange(N)=PGxyz(2); % Cross range from firing plane 
% Relative Velocity 
PV=[PV1; PV2; PV3]; 
PVR=VTR* aye (0; evs 10)" 
VDopp=PVR'? “PRxyz/SLR; *% Doppler velocity 
Vdop1(N)=VDopp; 
fd1(N)=2*VDopp/ lambda; % Doppler Frequency 
VRC: ,N)=PVR; 
VRr(N)=PV1; 


PVW=(PV1; PV2-WS; PVE-WE); 
HeadingR(: ,N)=VTR*( PVW/norm( PVW) ); 


a 


for i=l: N; 
VRMCi)j=norm( VRC: ,i)); 
VRrfCij=VRr(i)/VRMNCi); 
VRtCi)=C1-VRrfCi)})*VRMCi); 
beamD( i)=1*pi/180*SRR(i); 
end 
peaee 5 etime(clock,t2) 
ye cfr terior Gator br oe eee oe out jue data wrardstidedevedededededede ck Seeded sede rede sete ae ve 
diary ww 
Eee teen 
MN=MN 
ng=N 
MB=MB 
tg=t(N) 
ns eal 
dlatitude=(P2-G2)*180/pi 
dlongitude=PG3*180/pi 
Rg=GRange(N) 
FRg=FRange(N) 
CRe=CRange(N) 
CRemax=max( CRange) 
maxHeight=max( Height) 
maxfd=max( fd1) 
maxrange=SRR(N) 
VDOPmin=min(Vdop]1) 
VDOPmax=max( Vdop1) 
diary 
JO ETT wererereresevetease SAVE AND LOAD DATA FOR NEXT CALCULATION***7% 
PPE Ss OOmtr MAT 
SAVE 50010 MAXRANGE EXETIME] DR HEADINGR SRR HEIGHT T N FDI MB 
DEL AS50010 NAT 
SAVE 450010 TRAZ TREL VDOP1 VRM VRRF VRT VR VRR T SRR vei 


On vevetededs veisdevedeteeere i | —-Execute other program de vevevoveve decease Tesete te Wee eae a8 ee Tee 
7oPR 7o PRINT OUT TRAJECTORY 

*PRTRAK % TRACKING PARAMETERS & OUTPUT 

7440 7o RCS & RANGE ANALYSIS PROGRAM 

*RCSA *% RCS SUBFUNCTION 

*?RANGE3 % POWER REQUIREMENT 

“oI DATA 7 iden Matis, SURME ee 

OBE TTT TEE TE TE OEE OTE To TE ES EE TEE TE TEBE TE EEE ETE. TECTED TEE TEBE ETE PEPE PERE TEE EEE DEEP EEE BEER EER ER 
Ware e eee ee eee ete ee eee eee =r output of traj ectory*™ ieikikietetete eibte te teth ik tke ee 


'del GH. met 
plot(GRange,Height), grid; 
title('PROJECTIIE TRAJECTORY' ); 
xlabel(' ground range (m)'); 
ylabel(' projectile height (m)'); 
meta GH 

~Dresc 

%pause 
'del GC. met 
plot(GRange,CRange), grid; 
title('PROJECTILE TRAJECTORY’); 
xlabel(' ground range (m)'); 
ylabel('Cross range (m)'); 
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meta GC 
% pause 
Pap GES C 
!del RH. met 


plot(SRR,Height,'o'),grid; 
title('PROJECTILE TRAJECTORY’); 
xlabel( Slant Range (m)'); 
ylabel('projectile height (m)'); 


meta RH 
fe prtsc 
pause 
% !del TH. met 


paot(t,Height), 


xlabel(' flight 


grid; 
time (sec)'); 


ylabel('projectile height (m)'); 


meta TH 

% pause 

*pPrtsc 
!del TSGFC. met 
PlOtt(t ,oRange, 
xlabel( ‘flight 
ylabel('slant, 
meta TSGFC 

*pause 

presc 
!'del SGFCT. met 
plot(SRange,t, 
xlabel( ‘slant, 
ylabel(' flight 
meta SGFCT 

*pause 

omc Sc 
!del TVR. met 


t,GRange, t,FRange, t,CRange), grid; 
time (sec)'); 
ground, firing plane & cross ranges (m)'); 


GRange,t, FRange,t, CRange,t),grid; 
ground erorane plane & cross ranges (m)):; 
Gimemcsice a); 


femet(t,VRr, + ),grid; 


xlabel(' flight 


time (sec)'); 


ylabel( Vertical Speed'); 


meta TVR 
*pause 
Rpresc 

!'del TVH. met 


pilet(t,VRt,'+'),egrid; 
palabel( Flight time (sec) ); 
ylabel( ‘Horizontal Speed’); 


meta TVH 
*pause 
Presc 

!del FH. met 


plot(FRange,Height), grid; 
xlabel('ground range in firing plane (m)"); 
ylabel( ‘projectile height (m)'); 


meta FH 
*pause 
*prtsc 

!'del FC. met 


plot(FRange,CRange),grid; 
xlabel('Firing plane Range (m)'); 
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ylabel('Cross range (m)'); 


meta FC 
%pause 
PrEese 
!del GV. met 


plot(GRange, Adu eo ealvols 
xlabel( ‘Ground Range Gwe 
ylabel('Projectile Speed (m/sec)'); 
meta GV 
%pause 
Zp usc 
!del GRV. met 
plot(GRange, VRieb eer ide 
xlabel('Ground range (m)'); 
ylabel('Ratio of radial velocity to speed’); 
meta GRV 
ypause 
prt S¢€ 
Gor eee reser eK IK  TRAKING PARAMETERS CALCULATION AND OUTPUT * 2% 
yas Seid srisckie ak teris) se kor ete ears Tracking angular velocity wipe bicherive eke lttelerrel 
avell=diff(trel). /diff(t); 
avazl= oe ae i ialvesg sa ea)- 
Ake cob oe eta ok ih iy soedeve PLOT dedecdkdevededsdededsvedtdedededesdk deakaleedk as vdededk ve ves oe 
LOU SR 
Pal een); 
end 
del AVAEt. met 
plot(ta,avell, +. ).enid. 
title( ‘RADAR TRACKING ELEVATION ANGULAR SPEED '); 
Xlabel( In-flicht time tsce) 9) 
ylabel('Elevation Angular velocity (degree/sec)'); 
meta AVAEt 
% pause 
DEEse 
!'del AVAZt. met 
plot(ta,avazl, +'),grid; 
title('RADAR TRACKING AZIMUTHAL ANGULAR SPEED'); 
xlabel( In-=flicht timestsee) )- 
ylabel('Azimuthal angular speed (degree)'); 
meta AVAZt 
Ypause 


Owes ds iets sedvakieakdcdekiedkakacskkek wie Radar Tracking Parameter output dededeske desea aedk He seve 
'del tret. met 
PIOCGEVErei eri, 
title(' RADAR TRACKING ELEVATION ANGLE degree’); 
xlabel( ‘In-flight time (sec)'); 
ylabel( ‘Elevation Angle (degree)'); 
‘meta thel 
pause 
,apEeSe 
!'del trazl.met 
PlOBCRSEraz eye rns 
title RADAR TRACKING AZIMUTHAL PARAMETER’ ); 
xlabel('In-flight time (sec)'); 
ylabel( 'Azimuthal angle (degree)'); 


meta trazi 


%pause 


o/ 
/o 


% 
Zo 


prtse 

!'del AVAE. met 

plot(avell,'+'),grid; 

title(' RADAR TRACKING ELEVATION ANGULAR VELOCITY '); 
xlabel('MEASUREMENT BASE '); 

ylabel('Elevation Angular velocity (degree/sec)'); 
meta AVAE 

pause 

Brtsc 

!'del AVAZ. met 

plot(avazl,'+'),grid; 

title( ‘RADAR TRACKING AZIMUTHAL ANGULAR VELOCITY’ ); 
xlabel( 'MEASUREMENT BASE’); 

ylabel('AZIMUTH ANGULAR (degree)'); 

meta AVAZ 


*pause 
%prtsc 
EXETIMEl=etime(clock,t2) 


SRT TON TE TE TE TEE EET TE Te Ge Tee TE graph seacsesectakcctesescotvocscsesesesls sess seskectescseckestcsese st katt srsesk ok 


o 
ie) 


!del TVD. met 

Pitoe,t,fd]),egrid; 

title( DOPPLER FREQUENCY’ ); 
flabel{ Flight time (sec) ); 
ylabel('Doppler Frequency (Hz)') 
meta TVD 

pause 


Prtsc 


!'del treVD. met 

peor (tre! ,Vdopl),grid; 

title( DOPPLER Velocity’); 

xlabel('Radar tracking elevation angle (degree)'); 
ylabel( ‘Doppler Velocity(m)’') 

meta trehVD 


*pause 
ee esc 


!'del trazVD. met 

meot{ traz,Vdopl),grid; 

meitie( DOPPLER Velocity'); 

xlabel( ‘Radar Tracking azimuthal Angle (degree)'); 
ylabel( ‘Doppler Velocity(m)') 

meta trazVD 


*pause 
*Prtsc 


!'del RVD1.met 

mot ( ORK,Vdopl), grid; 

title( ‘DOPPLER Velocity'); 
xlabel('Slant Range(m)'); 
ylabel(' Doppler Velocity(m)' ) 
meta RVD1 


*pause 


ac sc 


!'del tVD1.met 
meot(t,Vdopl),grid; 
title({ DOPPLER Velocity’); 


xlabel('In-flight time(sec)') 
ylabel('Doppler Velocity(m)') 
meta tVD1 
PEESE 
Qote ae ese ce Here siedbocdescotsestieskocdkictiee Radar Tracking Parameter output Seavacsese desk dese feo a 
!'del tret. met 
plot(t, treljecrrds 
jeyie ces RADAR TRACKING ELEVATION ANGLE degree’ ); 
xlabel( In-flicht times sec rey. 
ylabel( ‘Elevation Angle (degree)'); 
meta thel 
% pause 
Dutse 
del trazl.met 
DilOtl te rraz ental. 
title( ‘RADAR TRACKING AZIMUTHAL PARAMETER’); 
xlapel( inofitene time see yay 
ylabel('Azimuthal angle (degree)'); 
meta trazl 


*pause 
Prcce 
ESET ETET EVE WORE Te Te VETS Sess Ie Te Te SE Ie TO Ie TE IE IE OT Te BIE TE IS Ie TEE IE UE OE Ie UE UE De IE 9S IE TE IE Ue 1K Ie IE IE Oh OC Oe Oe eae Ie Ie a ee 
Berens: RCS and Power requirement Estimation ** rrr 
load 50010 
to=¢ lock 
Bvl=1. 2*max(fdl); % video Bandwidth for lstutiiges 
Bv2=1. 2°2203]1;,; % video Bandwidth for 2nd filean 
es 1. 2812721; % Video Bandwidth for 3rd Jimiee. 

piace pee ee radar parameters calculation rok Lirh ara hae bk bt ae kd ead dea aie arise bore inak oro) or oe ey 
tases 100. 25: ye total loss 25405 
KT=4. O2Z1E-21; %o KI=K*T 
SN=10; % minimum signal to noise ratio 10dB 
Lamda=3E8/10525E6; % Wavelength meter 10525MHz 
Pn O02: “*o noise figure 2dB 
G=10"""4, 43; 7% TX Or RX antenna gain 44. 3dB 


corer ea "pi **3#KT*Fn*Loss*SN) /((G" bows, 
“ es! wlesieslectes'e os eee alton sles ao ate ve svs'ss ‘< ese + steels svexseveves evlesles'e este deste Jeveve se 


Zo Radar cross section & power racuipenens estimations 


for n=1:N 
if max(Height)==Height(n) % Summit point data 
ttop=t(n)- % time at summit 
Ltop=SRR(n); % range at summit 
end 
end 
Ute BURR REE EEE Switch J power requirement Wisds vee seat te seesay ce tiemage 


Bcon=Bvl*const; 
for n— lee 
SRRaC(n)=SRR(n); 
ta(n)=t(n); 
uD1=DRC(: ,n); 
h=HeadingR(: ,n); 
et=h'*uD1; 
rcS=rcsa(ct); + ECs mEnichlon 
eilCn j=eu: 
rcsl(n)=rcs; 


eee 


thetal(n)=acos(ct)*180/pi; 
L1=SRR(n); 
PtO=L1**4*Bcon/rcs; 


Pti(n)=Pt0; % power requied in free space 
if Li==Ltop % required power at summit 
ptop=Pt0; 

end 


Udore de dere se deve de de Fede ers Vege He HH Switene 2 power requirement Sesevedeslcsevedeseseteskeok vedo 
Bcon=Bv2*const; 
for n=1093: 1630 
SRRa(n)=SRR(n); 
ta(n)=t(n); 
uD1=DR(: ,n); 
h=HeadingR(: ,n); 
Se=n “udi; 
rcs=rcsa(ct); HoCeSuminc c 1On 
ctl(n)=ct; 
mes i(n)=rcs; 
thetal(n)=acos(ct)*180/pi; 
L1=SRR(n); 


PtO=L1****4*Bcon/rcs; 

Ptl(n)=Pt0; % power requied in free space 
if L1==Ltop % required power at summit 
ptop=Pt0; 

end 


Startech i=max( Ptl); 


end % end for loop 
8 Petes eTet ese ae Te 3rd filter SER llr rabensce d 
Bcon=Bv3*const; 
for n=1631: 1898 
SRRa(n)=SRR(n); 
ae 11 )=t(n); 
uD1=DR(: ,n); 
h=HeadingR(: ,n); 
ct=h'*uDl; 
rcs=rcsa(ct); meres Lunction 
@ci(n)=ct; 
resl(n)=rcs; 
thetal(n)=acos(ct)*180/pi; 
Li=SRR(n); 
PtO=L1**4*Bcon/rcs; 


been )=Pto; % power requied in free space 
if Ll==Ltop % required power at summit 
ptop=Pt0; 

end 


switch2=max(Ptl); 


end % end for loop 


y eee TORT RE REREE EE Out put Sar aerate DV ce Ta ee ce ROR tor earn atonal Gabon arrsn secs 
{del wwl 

diary wwl 

N=N 

flighttime=max(t) 


stoptime=max(ta) 
stoprange=max(SRRa) 
stopfd=fd1i( 1898) 
timetop=ttop 


rangetop=Ltop 
POWereep—pLtop 
EXETIME2=etime(clock,t3) 
TOTALTIME=EXETI MEW EXE PE 
diary 
Ore Meee Ieee TEE TE eae Tee Tee Te Te Hee He ae Ie graph Tadar CrOSS SECLCION WRK Rese weak sesekesksesk 
!del rcsll.met 
plot(SRRayresl. + eerace 
ylabel( ‘Radar cross section square meters'); 
xlabel('Slant range'); 
meta rcsll 
POLESC 
70 pause 
!'del rcsth. met 
plot(theralhres is — eric: 
ylabel('Radar cross section square meters ); 
xlabel('Aspect angle (degree) '); 
meta rcsth 
% pause 
Yoprtsc 
!'del rasp. met 
plot(ta,thetal, + ),grid; 
ylabel( ‘Radar Cross Section square meters s). 
xlabel( ‘Aspect angle (degree) ); 
meta rasp 
*o~pause 
Zpresc 
del rasl.met 
plottct lt oe) eri: 
Xlabel( incilaent timestcec)) 2. 
ylabel('Aspect angle (degree)'); 
meta rasl 
%eprtse 
oP ause 
9 vedtvedevods verre aie save Variable wewidkdedasedecedesises ode dea dese Heke desk de deve H ak se see see 
%: del Pw50010. mat 
“save Pw50010 ptimax ctl t fdl resi Li N thetal Peiesrna 
cia ripe tir Girone yea arante ro Ti ae eros prineoue power requirement 


ON OM EM CN ER ER CRIN IN ENR FL ELEN FV ER ODED EN 


oo ese 
'del powerl. met 
plot(SRRa,Pti,'+'),grid; 
title('RADAR POWER REQUIREMENT '); 
xlabel('Slant range (m) '); 
ylabel('Power requirement (Watts)'); 
meta powerl 


prtsc 

*pause 

wktkiskiekak RADAR CROSS SECTION CALCULATION Seve cdesedeslesesdevevenkeve a se vesevedk sc 

clear 

LOAD 50010 

for n=1:N 
WD1=DRie en) % radar to target unit vector 
h=HeadingR(: ,n); *unit heading vector of projectile 
ct=h'*uD1l; % antenna beam angle 


PEUNCt TOR mises. Cl] “RES Cec) 
% constant declaration 
b=0n a>, 2: % Shell base radius 
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l=Or5250; % shell cone(nuzzle) height 
i=0. 4144; % shell cylinder height 
ta=b/1; % =atan(b/1) i.e. half shell cone angle 
ta2=ta*ta; 
Gar—l/(l+ta2): 
ca=sqrt(ca2); 
Gx—adbs(ct ): 
@x2—CX CX} 
Sx2=(l-cx2); 
Sx=sqrt(sx2); 
resc=2*i*b*sx; 
b2=b*b; 
mes bD=pi*b2*cx; 
if cx>=ca 
mes=rcsctrcsb: 
else 
rt=ta*cx/sx; 
fee =rt rt; 
mese2—b2( asin( rt )trt*sqrt(l-rt2) )*cx; 
mes )=((b*1)/(sx2))*(1-cx2/ca2)***(3/2); 
fees=—roesit+rcs2+rcsct+rcsb/ 2; 
end % radar cross section case 2 
resO(n)=rcs; 
enn )=ct;: 
thetal(n)=acos(ct)*180/pi; 
end 
'del rcsll 
Meevot( SkR,rcs0, + ),grid; 
%ylabel('Radar cross section square meters’); 
*xlabel('Slant range’); 
“meta rcsll 
WPItsc 
“pause 
meet rcsth 
eect thetal,rcs0,'+'),grid; 
ylabel('Radar cross section square meters ); 
xlabel('Aspect angle (degree) '); 
meta rcsth 
*pause 
!'del rasp 
miler t,thetal,'+'),grid: 
ylabel(' Radar cross section square meters’ ); 
xlabel('Aspect angle (degree) ); 
meta rasp 
7opause 
!del rasl 
mer(ctl, + ),grid; 
xlabel('in-flight time (sec)'); 
ylabel('Aspect angle (degree)'); 
meta rasl 


PMeM eR Te Ke Ve HI Me He Fee end of program Kokovelvskoesdesdecrscscokscokaesdeaeske sks seok sek ot aesk ake sk ok kok 
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APPENDIX B. CROSS-SECTION ANALYSIS 


The radar cross-section has been discussed in Chapter 4. The derivation of the RCS 
equations is given below. The radar main beam direction is pointing in the direction 
D. The projectile axis (or heading) 1s pointing along h. Accordingly, the coordinate 
svstem (X,¥.Z) can be chosen so that x= D and the v-zZ plane is the plane on which the 
phvsical cross-sectional area of the projectile is to be projected. Another coordinate 
system (<1 2 cai be chosen omtliaus, — h and the projectile is a bodv of revolution 


around the x’ axis. These two coordinate systems can be set up jointly as follows : 


=D (B.1) 

i= D x h (B.2) 
|Dxh | 

f=afixD (B.3) 

pon (B.4) 


bh 
~~ 

I 

fir 

om 
S&S 
La 
- 


Tih Dir =D | 


=o ee (B.S) 
| xr | 

; Poe 

pe (B.9) 
ip 


Note that the coordinate (x.v.z) are used in this Appendix onlv and are not related 
to other parts of the thesis. 


The suface of the projectile is given by : 


ya! = f(x!) 


The RCS is the area enclosed by projecting of this surface on the y-z plane. 


ee CASE | 


(B.10) 


For | Deh|=>cosa, Figure 25 on page 38, the conoid 1s hidded in the shadow re- 


gion. Ihe projected areas can be easilv seen to be: 


Ay 


2 
. a 
3.14139 


0; | Dh | 


A, = 2ab 
ee 41 Deh | 
A, 
= 4,(D-h| 


(B.16) 


where a 1s the diameter of the projectile and bis the cvlinder height. The total projection 


area in: 
Te 0, + 0, + 6; 


Pee CASE 2 


(B.17) 


When |D-A|<coso, the geometry of the target 1s shown in Figure 26 on page 


Be liie RCS value consists of four different parts: A, , 4,. A,. A, . 


below : 


For area 4;: 


' hee. 2 
= (x’—x,') tan Oe 
where 


c = tan o 
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Each of them ts given 


(B.18) 


(B.19) 


(B.20) 


2 2 2 2 2 
(C Slane gen) tok, eae ae eee tan’é 


; 
? as 
Xq (anauiee 
(Cae ee aaa) 


ee. 
| es aN ee 


ce ae 5 
(lance Cy) 


where x,’, ¥, 1s the maximum point of the ellipse. 


Z 
C 


tana 


a] 
4 
ba ~ 


OG us ea ane eee 


At a pom which 1s pp, x = 15 then: 


Nene epee a 


or 
; 
eae eae 
(ames, 
a — Le’? 
SO the coordindie at 7,45 (-- = 
Lata ae 

Ae) 

CE Ae. pe 
Se ES 
MENON Care 
or 
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ayn = 


~O ) On the aac 


(B.28) 


(B.31) 


oe 
A, = 2L(——> is (B.32) 
(ele es 
2 2 
ane — Cc) 
4 asp 2 
= L*c( tan*Oo6, +¢ re ee (B.33) 
Laneere 
mete, be the normial unit vector of surface 4; 
A 4 “A 
re OCS elie 
SS (B.34) 
Pe 2 4 me 2 

Meo) 07, COS Oo + Sin & COs U,., 

then the projection area for 4A, 1s: 
au 4 19 2 D, pees A 

L* tan o( tan°6,,, + tan a)'!"( tan’6,,, — tan’a)'/(D +N) 

a an ae 7a (B.35) 
eI) Boge 
cos"@ 
a esp 375 

= 7 - = (B.36) 

SUED sy EG sno 


A, is part of a semi-circle (Figure 26 on page 39). A straight forward intergration gives: 


2, 2:7!) la tan % tan-9 
a —ai{isn (—— )+—— /!1-—} (B.37) 
: | tan ep | | tan oe | lor Oe 


ee! tan 9 tan % tan-o a 
rm) ee (B.38 
: ( aT Fe Gane. | tan’ 65, | , 
A, 1S a rectangle: 
A, = 2ab UB ee 
SO 
A, | D xh | (B.40) 
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where a is the radius of the projectile base, and 6 1s the height of the cylindrical part of 
projectile. 


A, 1S a semicircle, hence : 


A, 


2 
a 
3.14159 — (B.41) 


Aan (B.42) 


I 


O4 
The total cross-section ¢, 1s the sum of upper 4 projection areas. 


0, = 6,+6,+6,+0; (B.43) 
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APPENDIX C. PROPAGATION FACTORS 
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. 0000 
. 6000 
. 4000 
: 3000 
. 3000 
. 4000 
. 1000 
. 4000 
. 7000 
. 2000 
7000 
2 000 
. 7000 
. 0000 


. 6000 
. 8000 
. 9000 
. 6000 
. 3009 
. S008 
. LOGE 
. 5000 
. 0000 
. /000 
. LOGE 
. 1000 
. 8000 
, 2060 


APPENDIX D. MISSED DETECTION AND FALSE ALARNI 


clear 

MN=1000; 

Npeak=MN/ 2; 

dtmax=0. 08650; 

Bv 1=63900; 
FNnoise=Bvi*dtmax/Npeak; 
SNRO=10. 0; 

seep— 0.5; 

Nf=5; 

Mf=29; 

thshstep=0. 005; 

thsh0=0. 355; 

for m=1:Mf 
thsh=m*thshstept+thsh0; 
thresh(m)=thsh; 
FNEnoise(m)=FNnoise*2*asin(sqrt(thsh) )/pi; 
end 

X=SNRO-step; 

for n=1:Nf 

x=xtstep; 

SNR(n)=x; 
segnr=10**(0. 1*x); 

% transfer the signal-to-noise parameter to be used in fdetect.m 
snr=sgnr 

'del hm tempsnr. mat 
save hm tempsnr snr 


for m=12 Mt 
thsh=thresh(m); 
pfan=exp(-thsh*sgnr); 
Pfa(n,m)=pfan; 
FNerror(n,m)=pfan*FNEnoise(m); 
end 
Pmiss(n,1)=Pmissdct(sgnr,0,thresh(1)); 
for m=2:Mf 
Pmissd=Pmissdct(sgnr,thresh(m-1),thresh(m) ); 
Pmiss(n,m)=Pmiss(n,m-1)+Pmissd; 
end 

end 


!'del hm pfaulse. mat 
save hm pfaulse Nf Mf thresh FNEnoise SNR Pmiss FNerror Pfa 


'del hm pfaulse 
diary hm pfaulse 


disp(' Threshold Missed Noise Faulse' ) 
Grsp( Detection Error Alarm = > 
for n=1:Nf 


SN_ratio=SNR(n) 
(thresh; Pmiss(n,:); FNerror(n,:); Pfa(n,:))'; 


~3 
{sy 


end 
diary 
pause 


!'del hm pfaulse. met 


plot (thresh, Pmiss(1,:), + , thresh, FNerrorq@ie), © mecsiGe 
title( 'Signal-to-Noise Ratio = 10.0 db’) 
xlabel('Threshold-to-Signal Power Ratio’); 

ylabel('Missed Detection (+) & Noise Error (0)') 

meta hm pfaulse 

pause 


plot (thresh, Pmiss(2,:), + , thresh, FPNerror(2..)) © omeaice 
title('Signal-to-Noise Ratio = 10.5 db’) 

xlabel( 'Threshold-to-Signal Power Ratio’); 

ylabel( ‘Missed Detection (+) & Noise Error (0)') 

meta hm pfaulse 

pause 


plot (thresh, Pmiss(3,:), + , thresh, FPNerror(3.9 = © jmeoiee 
title( Signal-=to-Noise Ratio = tle 0db.) 
xlabel('Threshold-to-Signal Power Ratio’); 

ylabel( ‘Missed Detection (+) & Noise Error (0)') 

meta hm pfaulse 

pause 


plot (thresh, Pmiss(4,:), +’, thresh, FNerror( 4): jo seen 
title( Signal-to-Noise Ratio = 11.5 db } 
xlabel('Threshold-to-Signal Power Ratio’); 

ylabel( ‘Missed Detection (+) & Noise Error (o0)') 

meta hm pfaulse 

pause 


plot (thresh, Pmiss(5,:),'+', thresh, FNerror(5, 3) omemuen 
title( Signal-to-Noise Ratio = 12) 0edp 

xlabel( 'Threshold-to-Signal Power Ratio ); 

ylabel( ‘Missed Detection (+) & Noise Error (0)') 

meta hm pfaulse 

pause 


Vertovedevededtedetedevese sete vee Seve tect se cece Fe Fe eer e Se De eee Teen TOTO TOE Te TET TE Ue Oe TOE TONE TC Te cee ee er en 
function pmsd = Pmissdct(snr, thsha, thshb) 

% The parameter snr for the singal-to-noise ratio (NOT in db) should 
% have been saved to hm tempsnr. mat. 

%tdel hm tempsnr mat 

*%save hm tempsnr snr 

pmsd=quad(' fdetect’ ,thsha,thshb,0. 0001); 


aloutewonte aloaten’ouiouts ana anlon nace evans woaleneelculcatswonontewosioweniseisn'celenianiowenoeweaieweutsuts Fo gts Joute low'e tects Se ctes’e de Je Jove fo wie cle fe ode cto ge ses'eclos'c 
Gt 4h £O FL OR CS £H FRX £% BR OB FE FX CR EX FX FE FR FY FR CR CX CHT FT EX FL FR 4S FR ET ER FLFR FR FR eeeuneuevenen as NER ER ON RENAN EN EN AN CN EN EN EN OW ON GR ENED ee s . 


ae ou ev ee 


% The parameter sgnr has to be specified in the program. 
function pdens = fdetect(v) 

load hm tempsnr 

ZAV=2 Sis Ley aa, 

LIloz=log( BESSELN(0,ziv)); 
Lpdens=LIiIoztlog(snr)-snr*(v7+1); 
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pdens=exp( Lpdens ); 
Wwrevesevedeseveskiewe keke end of program dere de dese Teale He de ae ae ae He ae ae de ae ae ve ae oe ae Teak ae ae ak Te ake ae Me Ie He Ie 
Uprededeve de vevore deve Hake ve He se He output data wkecadededevesdtdcdedvabsededevedevedese deve ve secede ae Fee He ae vere Hee se ve 
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